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ABSTRACT 

In a study of the planetary-scale vorticity budget of 
the northern summer monsoon, Holton and Colton (1972) found 
that a strong vorticity sink is needed to balance the gener- 
ation by horizontal divergence at 200 mb. In this work we 
use numerical experiments to examine the hypothesis that 
the fluctuation of the divergence forcing as observed by 
several studies serves as this sink. A 3-level numerical 
model is used with the thermal forcing specified by the hori- 
zontal divergence distribution observed by Krishnamurti 
(1971). Two experiments are carried out, one with heating 
steady in time and the other fluctuating in time with a 
period of 10 days. The time-mean fields of the quasi-steady 
state solution in both experiments show a significant west- 
ward phase shift of the Tibetan high, thus the fluctuation 
of heating cannot be a sufficient sink. On the other hand, 
the positions of the simulated tropical upper tropospheric 
troughs agree quite well with those observed. This sug- 
gests that damping due to cumulus transport is more likely 
to account for the vorticity sink. Transient wave activi- 
ties are also found at the upper level in both experiments. 
They appear to be drawing energy from the planetary-scale 
flow and could serve as a vorticity sink. These waves are 
most active in the fluctuating heating case, suggesting 
that local barotropic instability at upper levels may be 
more important than that implied by the time-mean wind. 
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I. INTRODUCTION 



In the study of the northern summer planetary-scale mon- 
soon circulation, Holton and Colton (1972) discovered that 
very strong damping is necessary to balance the budget of 
time-mean planetary-scale vorticity at the upper levels. 
Specifically they found that without this strong damping to 
counterbalance the generation of vorticity by horizontal 
divergence the amplitude of the long waves of the monsoon 
circulation would be too large at the 200 mb level and the 
phase would be shifted westward by up to 5000 km. 

This same phase problem was also apparent in a numerical 
study by Abbot (1971), who used a 2 layer quasi -geostrophi c 
model with prescribed heating to study the monsoon. 

In looking for a source of this strong damping mechanism, 
Holton and Colton suggested two hypotheses which may explain 
this strong damping. The first is the vertical mixing of 
vorticity and momentum by strong cumulus updrafts and down- 
drafts. The second hypothesis is related to the observed 
time fluctuation of the monsoon system. If the large-scale 
divergence contains a strong fluctuation then a correlation 
term of the divergence and vorticity fluctuations will ap- 
pear in the time-mean vorticity equation and may reduce the 
generation of the time-mean vorticity by the time-mean 
divergence. 
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Murakami (1974) used a three-component, quas i -s pectral 
model with prescribed heating which has a fluctuating compo- 
nent. He found that it can reasonably simulate the ampli- 
tude and the phase of the long waves associated with the 
northern summer monsoon circulation. But his model suffers 
from the lack of the zonal basic flow and the resultant ver- 
tical profile of the divergence field differs considerably 
from observations. Therefore it is not clear which mecha- 
nism causes the seemingly correct phase and amplitude of 
his model . 

The purpose of this study is to examine the second hypo- 
thesis of Holton and Colton using the non-linear, finite 
difference model of Monaco and Williams (1975). This same 
model was used to simulate the northern summer monsoon cir- 
culation by Beilis (1975). In his study, Beilis found that 
the prescribed heating, which was specified by the 200 mb 
mean horizontal divergence field observed by Krishnamurti 
and Rodgers (1970), can reasonably excite several important 
features of the time-mean monsoon circulation. However, if 
extremely large values of damping are not included, a phase 
shift problem still exists, although the amplitude of the 
disturbance seemed to agree with that observed. 

In this study we will carry out two controlled experi- 
ments. One will have the heating specified to be independ- 
ent of time and the other will have a fluctuating heating 
component similar to that observed and used by Murakami. 
Hopefully by comparing these experiments we will be able to 
evaluate the effects of this fluctuation in the monsoon heating. 
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Another purpose of this study is to simulate the baro- 
tropic instability of the tropical easterly jet associated 
with the northern summer monsoon. This instability arises 
as a result of the strong easterly jet in the upper tropo- 
sphere and has been studied by Colton (1974) and Beilis 
(1975). It represents an important seal e- i nteracti on mecha- 
nism for the monsoon circulation, according to the observa- 
tional study by Krishnamurti (1971). 



II. THE NUMERICAL MODEL 



A. BASIC MODEL 

The numerical model used was formulated by Monaco and 
Williams (1975) and is similar to that described by Arakawa 
and Mintz (1974) who used the primitive equations in spheri- 
cal coordinates with "sigma" as the vertical coordinate. 

In this study the domain of calculation is restricted to 
the band between 18°S and 46°N for the sake of computer 
economy, and a free slip condition is applied at the lateral 
boundaries. The model has been previously used by Beilis 
(1975) for the area of 0°E to 180°E only. 

B. THE GRID SYSTEM 

The horizontal distribution of variables is illustrated 
in Figure 1 which is the so-called C-scheme arrangement 
(Arakawa and Mintz, 1974). The grid spacing is 4 degrees in 
both latitude and longitude. The troposphere is divided 
into seven equally spaced sigma layers in the vertical with 
three main reporting levels, as shown in Figure 2. The 
sigma coordinate is defined as 




where p is the pressure of the sigma level, and p t is 
the constant pressure height of the tropopause set at 100 
mb, and tt is the terrain pressure defined as 

tt = p s - P t (2.2) 
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The surface pressure, p $ , was set at 1013 mb initially. 

The vertical boundary conditions are ^ = 0 at a = 0, (p 
= p t ) and a = 1 , (p = p $ ) . 

C. SPECIFICATIONS OF HEATING AND DAMPING 

To evaluate the effect of fluctuating heating as sug- 
gested by Holton and Colton (1972), two controlled experi- 
ments are carried out. In both experiments the heating func- 
tion, H, contains two components: 



H 



H s + 



'AS 



(2.3) 



where H $ , the zonally symmetric heating, is given by a 

continuous adjustment to a specified equilibrium temperature, 
★ 

T , using the following formula. 



H 



s 





(y »q )j 



(2.4) 



Here T is the temperature and x is a specified adjust- 
ment time set to two days in this study. The equilibirum 
★ 

temperature, T , is prescribed to maintain the zonal mean 
north-south temperature gradient. 

The horizontal distribution of heating is derived from 
the 200 mb time-mean divergence field observed by Krishna- 
murti (1970) as follows: 

Q(°K day"^) = Divergence (sec~^) x 10^ (2.5) 

The field is smoothed on the boundaries by reducing the 
strength of the heating by 50% at the points just inside 
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the boundaries and 80% at the boundary. The entire field 
was smoothed by a five point space smoother in order to re- 
move small scale features. In the vertical, the full value 
of heating derived by (2.5) is specified at level 2 and a 
vertical distribution is assumed such that the heating at 
level 1 and 3 are 50% and 25%, respectively, of the value 
at level 2. 

The zonally asymmetric component, takes the follow- 

ing form, 

H a$ = Q x (1 .0 + A ( S I N y- t)) (2.6) 

where A is the amplitude of the heating fluctuation whose 
period is T. In the first of the two experiments A = 0; in 
the second, A = 0.5. A 10 to 15 day period was observed by 
Murakami (1974) and Krishnamurti and Bhlame (1976) in the 
northern summer monsoon. Here ten days is chosen for the 
periodicity T for the sake of computation economy. 

Linear damping terms are included in the momentum equa- 
tions to represent the effect of scale- independent friction. 
Horizontal dif fusion is also introduced to include the effect 
of scale-dependent friction for it damps the small scale 
fields, which are generated by the finite difference scheme 
in the model . 

D. INTEGRATION SCHEME AND INITIAL CONDITIONS 

The time integration is comprised of continuous sections 
of one Matsuno step and four leapfrog steps, with a time in- 
crement of six minutes. 
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At the onset the pressure levels were initialized at 
standard heights and temperatures were also specified by 
the zonally averaged climatology. The model was started 
without the diabatic forcing function until a quasi-steady 
circulation was established. At this point the diabatic 
heating was introduced. For both experiments the model 
runs for 25 more days from the quasi-steady state with the 
data of the last 10 days used for computing the time-mean 
fields. All data presented is interpolated to constant 
pressure surfaces of 250 mb, 500 mb, and 850 mb. 
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III. RESULTS 



The principal circulation features of the northern summer 
monsoon are the Tibetan high, the Pacific tropical upper 
tropospheric trough (TUTT), the tropical easterly jet, the 
Mexican high and the Atlantic TUTT. The simulation of these 
features by the model are compared with observations, mainly 
those by Krishnamurti and Rodgers (1970). 

A. TIME MEAN FIELD 

1 . Steady Heating 

The 10 day time -me an fields of wind, geo potential 
height, temperature and vertical velocity are shown in Fig- 
ures 4-7. In Figure 4a the 250 mb wind field shows an anti- 
cyclonic wind center near 30°N, 52°E which resembles the 
observed Tibetan high although the position of the generated 
feature is shifted 30° to the west. The calculated Pacific 
TUTT extends from 120°W, 40°N to 180°W, 20°N, approximately 
the same position as the observed TUTT. The TUTT has a 
northeast-southwest tilt and the amplitude is slightly greater 
than that observed. The strongest easterly jet is developed 
by the model in the southern periphery of the Tibetan anti- 
cyclone. The jet maximum is located at approximately 60°- 
90°E, 10°N, with a speed of 22 ms ^ which is an underestimate 
of the observed maximum speed of ~30 ms - ^. 

The Mexican anticyclone is centered near 130°W, 25°N. 
This feature is again shifted to the west of the observed 
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position by 25°. The amplitude of the disturbance compares 
well with the observed. 

The Atlantic TUTT, situated in the mid-Atlantic in 
the area of 15°N-40°N, 20°-60°W is well developed but again 
the amplitude is slightly larger than that observed. As for 
the Pacific TUTT, there is no apparent phase shift problem 
when compared to observations. 

There are no good observational data with which to 
compare the 500 mb flow shown in Figure 4b. The 850 mb flow 
(Figure 4c) clearly shows the low-level cyclonic flow below 
the upper-level Tibetan and Mexican anticyclones. Westward 
phase shifts of circulation centers also occur at low level 
and is consistent with the upper level flow. 

In Figure 5a, the 250 mb geopotential, the Tibetan 
and Mexican highs are the main features. The Pacific and 
Atlantic TUTT s are weak in comparison to observations but 
well defined. 

Figure 6a, the 250 mb temperature field, shows that 
the Tibetan and Mexican highs are warm core systems. The 
TUTTs in both the Pacific and Atlantic are clearly indi- 
cated by cold regions. These temperature anomalies are 
present throughout the vertical in both the 500 mb field 
(Figure 6b) and the 850 mb field (Figure 6c). 

As expected, Figure 7, which shows the 400 mb verti- 
cal velocity closely resembles the horizontal diabatic heat- 
ing distribution. This is consistent with the warm core 
energetics of the monsoon in which the long wave kinetic 
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energy is converted from the available potential energy 
generated by diabatic heating. 

Figures 8 and 9 show the mean divergence and vortic- 
ity fields, respectively. The 250 mb divergence (Figure 8a) 
is a reasonable reproduction of the 200 mb divergence field 
as observed by Krishnamurti . This indicates that the 
specified heating function is properly chosen for simulating 
the forcing. The vorticity fields are representative of 
the wind fields, shown in Figure 4. A negative vorticity belt 
surrounds the globe from 8°-25°N. The two TUTTs are clearly 
shown by fingers of positive vorticity extending into the 
negative vorticity belts. 

2 . Fluctuating Heating 

The 10-day mean fields of the fluctuating heating 
forcing function are Figures 10-15. Only the 250 mb level 
of all the fields is shown except in the cases of vertical 
velocity (Figure 13) and the time-mean wind fields (Figures 
10a, b,c). A close comparison of these fields with the time- 
mean fields for the steady heating case could not reveal any 
significant differences between the two cases. The phase 
shifts of the Tibetan and Mexican highs are almost exactly 
duplicated by the fluctuating heating, resulting in almost 
identical amplitudes and phases as those for the steady 
heati ng . 

The above results suggests that the effect of 
fluctuating divergence forcing in the summer monsoon suggested 
by Holton and Colton is not prominent enough to explain the 
phase problem in the diagnostic simulation of the time-mean 
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field. This conclusion is true if the amplitude of the fluc- 
tuating component in our model is representative of the 
actual observed amplitude. Since limited observational know- 
ledge (based on Kri shnamurti , 1970) does indicate that the 
actual amplitude of fluctuation may not be much larger than 
what is specified here, Holton and Colton's hypothesis seems 
unlikely to be valid. 

B. TRANSIENT FIELD 

The second purpose of this study is to examine the synop- 
tic-scale transient disturbances associated with the plane- 
tary-scale monsoon circulation. 

1 . Steady Heating 

Figure 16 is a time-longitude section for the 12-hour 
interval 200 mb vorticity field at 12°N from day 15 to day 25. 
In the western hemisphere from 0°-90°W there is evidence for 
strong synoptic-scale westward-propagating wave activity. 

The wavelengths of these disturbances vary from 3500 km at 
30°W to 2200 km at 60°W. In some cases the wavelength can 
be as short as 1000 km at 180°W to 90°E where the disturbances 
appear to be much more irregular with long wavelengths at the 
beginning of the time series and sporatic, short wavelengths 
toward the end. 

Upstream from the jet maximum near 50°E the disturb- 
ance activity is very minimal with small amplitude. The waves 
begin to grow in the area of the easterly jet maximum and 
continue to grow as they propagate westward. Downstream of 
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the easterly jet maximum the synoptic waves appear to be 
much more organized comparing to upstream of the jet. 

Krishnamurti (1971) and Colton (1972) have found 
that the upper-level synoptic waves are largely due to the 
barotropic instability associated with the easterly jet. 

To examine this possibility as the generation mechanism of 
the waves, observed in this model, we plot the meridional 
gradient of the basic absolute vorticity 



where 3 is the gradient of the coriolis parameter and u 
is the time-mean zonal wind component, as a function of x 
and y in Figure 17. It is obvious that the sign of this 
quantity changes within the north-south channel throughout 
the domain. Thus it is quite probable that the barotropic 
instability exists everywhere at the 250-mb level. The im- 
portant role played by the easterly jet in this regard is 
reflected by the fact that the entire latitude band of 12°N 
has negative values for this quantity. 



stream of the easterly jet and increases to a maximum down- 
stream further suggests that the mean kinetic energy is con- 
verted into eddy kinetic energy in the vicinity of the 
easterly jet maximum. There should be a stable area upstream 
of the easterly jet maximum where the amplitude of the dis- 
turbances is small. The amplitude near the jet maximum may 
be small but the change in amplitude is large as the wave 




2 - 



The result that the wave amplitude is minimum up- 
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propagates past the jet indicating that the jet maximum is 
the most unstable region. This result indicates that the 
most active area for synoptic-scale wave disturbances is in 
the Atlantic, Caribbean and Central America, not in the 
Pacific and Indian Ocean areas. There are two possible 
factors in the model that influence this: (1) The model 

position of the Tibetan high is shifted by 30° to the west 
from the observed position, and there is a corresponding 
shift in the jet maximum from the Indian to the Arabian Sea. 
The most active wave region should be downstream from the 
jet maximum so this region is shifted from the Arabian Sea 
and Africa to the Atlantic area. Also the Atlantic TUTT 
has been known to be an active area for upper-level cold core 
lows to develop. It is possible that the synoptic waves 
observed here are also a reflection of the cold core lows. 

On the other hand, the absence of strong disturbance acti- 
vity in the Pacific TUTT region, in contrast to the finding 
by Colton (1972), may also be explained by the westward 
phase of the easterly jet maximum. (2) This model contains 
a horizontal diffusion which tends to reduce the zonal wind 
shear and in so doing the barotropic instability. This dif- 
fusion was not used by Colton (1972). It is noteworthy, 
however, that even in this model the barotropic instability 
mechanism is still found to be prominent downstream from the 
jet maximum. An optimum use of the diffusion, to alleviate 
the problem associated with the finite differencing and at 
the same time avoid strong artificial smoothing of the basic 
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wind shear, should give us a better simulation of this 
mechanism. 

As an example for the barotropi ca 1 1 y unstable wave, 
a sequence of 12-hour interval 250 mb wind analysis is shown 
in Figures 18 a-d, from 0000 GMT, day 22 to 1200 GMT day 23. 
Here we shall follow the development of three waves in time. 
Their positions at 80°W, 4°-16°N; 60°W, 4°-16°N and 30°W, 
4°-16°N on 0000 GMT day 22 are indicated by dashed lines. 

As time progresses each wave propagates to the west with a 
wavelength of -3500-4000 km. The meridional tilt of the de- 
pi cted waves i s approximately northwest-southeast. From 
Figure 4a, the maximum time-mean easterly zonal wind is near 
12°N in the vicinity of the waves in question. Thus in the 
12°-16°N band the time-mean zonal wind has a northeast-south- 
west shear which opposes the wave tilt. This is clearly 
another indication that the waves are ba rotropi ca 1 1 y un- 
stabl e . 

2 . Fluctuating Heating 

Figure 20 shows the time-mean meridional gradient of 
the absolute vorticity for the fluctuating heating case and, 
as expected, strongly resembles Figure 17 for the steady heat- 
ing case. There is no question on the plausibility for baro- 
tropic instability to exist. 

The time-section of the mean vorticity. Figure 19, 
also points out that the greatest wave activity is downwind 
of the maximum jet. The amplitude of the waves, however, 
appears to be greater in this fluctuating heating case, though 



23 



the range of wavelengths remains the same as the steady heat- 
ing case. There is also no noticeable change in phase speed 
of ~7.6° day”^. This increase in amplitude is readily 
apparent in the waves generated in the Atlantic area of 20°- 
90°W, 8°-16°N. This result indicates that the fluctuation 
of the total forcing field may actually enhance the barotropic 
energy conversion. In Figure 21, by following the waves 
initially located at 85°W, 8°-16°N; 55°W, 10°-16°N and 30°W, 
8°-16°N; it is found that the tilt of the waves, northwest- 
southeast, again opposes the horizontal shear of the basic 
wind. 
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IV. CONCLUSIONS 



One major conclusion of this study is that the westward 
phase shift problem of the simulated Tibetan and Mexican 
highs in the absence of very strong damping is not likely to 
be explainable by the influence of the fluctuating heating. 

On the other hand, the tropical upper tropospheric troughs 
do not exhibit this phase problem. This may indicate the 
importance of damping due to cumulus transport of momentum 
and vorticity because cumulus convection is predominant below 
upper level high pressure areas only. 

In this model the upper-level barotropic instability 
seems to be quite important for most of the area, even though 
the easterly jet was not developed to the fully observed 
strength, perhaps due to the inclusion of the somewhat 
strong horizontal diffusion. This instability seems par- 
ticularly important when time-fluctuations of the planetary- 
scale is included through the imposed fluctuating heating. 
Thus the local barotropic instability at upper levels may 
be more important than that implied by the horizontal shear 
of the time -mean wind. Further investigations on this 
mechanism would be desirable to see whether the observed 
upper level disturbances, such as the cold core lows, are of 
this nature. Incorporation of cumulus parameterization and 
topography in the model should further enable us to study 
the possible implication of this upper level barotropic in- 
stability in the vertical development of other tropical 
disturbances. 
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In addition to the phase shift problem, Holton and 
Colton (1972) also found very large amplitudes for long 
waves in the case of weak damping in their linear model. 
Our results do not have this amplitude discrepancy. This 
indicates that the transient waves due to non-linear baro- 
tropic energy conversion serves as an energy sink of the 
planetary-scale flow. 
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Interval is 1 mb hr- 
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Figure 16. 250 mb vorticity time-longitude section at 

12°N, steady heating. Hatched area indicates 
positive values. 
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Figure 19. 250 mb vorticity time-longitude section at 

12°N, fluctuating heating. Hatched area 
indicates positive values. 
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fluctuating heating. Hatched area indicates 
positive values. 
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